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ABSTRACT 



We report mid-infrared (MIR) imaging observations of the carbon star U 
Ant made with the Infrared Camera (IRC) on board AKARI. Subtraction of the 
artifacts and extended PSF of the central star reveals the detached dust shell 
around the carbon star at MIR wavelengths (15 and 24 /xm) for the first time. 
The observed radial brightness profiles of the MIR emission are well explained 
by two shells at 43'' and 50'' from the central star detected in optical scattered 
light observations. Combining Herschel/P ACS, AKARI/FIS, and ^i^^i?//IRC 
data, we obtain the infrared spectral energy distribution (SED) of the thermal 
emission from the detached shell of U Ant in a wide infrared spectral range of 15- 
160 /im. Thermal emission of amorphous carbon grains with a single temperature 
cannot account for the observed SED from 15 to 160 /im: it underestimates the 
emission at 15 /im. Alternatively, the observed SED is fitted by the model that 
amorphous carbon grains in the two shells have different temperatures of 60 and 
104 K, which allocates most dust mass in the shell at 50". This supports previous 
suggestion that the 43" shell is gas-rich and the 50" one is dust-rich. We suggest 
a possibility that the segregation of the gas and dust resulting from the drift 
motion of submicron-sized dust grains relative to the gas and that the hot dust 
component associated with the gas-rich shell is composed of very small grains 
that are strongly coupled with the gas. 

Subject headings: stars: AGB and post-AGB — stars: carbon — stars: mass-loss — 
infrared: stars 
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1. Introduction 

Mass loss from Asymptotic Giant Branch (AGB) stars is one of the major steUar 
sources that supply dust grains to t he interstella r medium (ISM) and is a key process of 



the material circulation in galaxies ( jPraine 



sources even in low metallicity environments ( ISloan et al. 



20091). AGB stars a re suggested to be efficient 



20091 ) . Mass loss is also crucial 



to the late evolutionary stages of low- to intermediate-mass stars. Radiation pressure on 
forming dust grains is thought to play a significant role as a mass-loss mechanism, but there 
still remain several open issues, such as the trigger of the mass-loss process, the dependence 
of the mass loss on the stellar parameters, and the change of the mass loss along the 



stellar evolution ( IHabing 



19961 ). Far-infrared (FIR) studies based on /i?74S' observations 



indicate the presence of detached shells around carbon stars, suggesting that the mass-loss 



process is sporadic rather than stea dy due to thermal pulses in the star ( jWillems fc de Jong 



19881 : 



van der Veen fc Habing 



19881 ). Detailed an alyses of FIR data in fact confirm the 



presence of detached shells in several AGB stars (Waters et al.. .1994 : .Izumiura et al. 



201 



Ueta et al. 



20101 : 



Ladial et al 



in detail with CO line observations ([Yamamura et al. 



Linda vist et al. 



scattered light ( 



1999 



Schoier et al 



19961 . 



201 oh. Detached shells are also invest i gatec 



Gonzalez Delgado et al 



19931 : 



Qlofsson et al. 



19961 . 



20001 : 



2005h and their 



2m 



200d; 



spatial structures are studied in optical 



Maercker et al 



2010h . While CO and 



optical scattered-light observations provide detailed information on the shell dynamics and 
structure, infrared observations yield estimates of the dust mass and temperature of the 
detached shell. Mid-infrared (MIR) data can add crucial constraints on the temperature 
and mass of dust grains because they measure the Wien side of the thermal emission from 
the detached shell. However, no detached shells of AGB stars have so far been detected in 
the MIR. 



U Ant is an N-type carbon star of irregular variability at a distance of 260 pc 
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(KnaDD et al. 


2003 


). 


(Izumiura et al. 


1997) 



nvelope around U Ant is studied in the 



19971 ). th e CO emission (|01ofsson et al. 



the optical scattered light f iGonzalez Delgado et al. 



2001 



1996; 



2003 



Maercker et al. 



Maercker et al 



IR 

20101), and 

X 

2010). These 



studies resolve 5 shells at - 25^^ 37^^ 43^^ 50 



2, 3, 4, and 5, respectively), 
for the existence of full she^ 



1^^ . and 3^ from the central star (hereafter shell 1, 



Maercker et al. 



( 120 lOl ) argued that there is no strong evidence 



with PACS (IPoglitch et al. 



s insi de shell 3. Very recently the shell aroun d U Ant is stu died 



20101) onboard Herschel Space 



2010), which resolves the detached shell clearly in the FIR (IKerschbaum et al 



Observatory ( 



Pilbratt et al. 



20101). In this 



Letter^ we report the first detecti on of the detached shell around U Ant at MIR waveleng ths 



with the Infrared Camera (IRC; 



Onaka et al. 



20071 ) onboard AKARI flMurakami et al. 



20071 ). The observations are made under the AKARI mi^^ion program "Excavating Mass 



Loss History in Extended Dust Shells of Evolved Stars" (MLHES, PI: I. Y.). 



2. Observations and Data Reduction 

Observations of U Ant were made in the L15 (15 /xm) and L24 (24 /xm) bands of the 
IRC on 2006 December 20. The FWHM of the point spread function (PSF) is 4'.7 and 6'.'8 
for L15 and L24, respectively. The obtained dataset consists of long- and short-exposure 
images. Since U Ant is very bright in the MIR (^ 100 Jy at 15 /xm), the long-exposure 
images are all saturated around the central star position and severa l artifacts appear in the 
images due to refiections between the optical elements of the IRC ( ILorente et al.ll2008l ). 
We use the short-exposure data to complement the saturated inner region and employ the 
long-exposure data to investigate the outer region. The short-exposure data of L24 are not 
saturated even at the central star, but those of L15 are nearly at the saturation level. Thus 
the very central part (r < 7") of the L15 image is not used in the following analysis. The 
short- and long-exposure data are smoothly connected at the unsaturated region by scaling 
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with the exposure time, (a factor of 28). 

Accurate determination of the artifacts and the PSFs with a large dynamic range are 
needed to detect faint emission from the detached sheU. We employed the data of the bright 
star IRAS 22396—4708 taken in the MLHES program as a calibration source, the archival 
data of CRL618 of an Open Time program, and IRAS F06009— 6636 taken in the IRC 
calibration observations to investigate the artifacts and PSFs. The first two objects were 
used to investigate the artifacts and the PSFs in the outer part, while the last one was used 
to determine the inner PSFs, where the former data were saturated. The artifacts show 
movement patterns different from real objects in the image. Using this characteristic, we 
separated the artifacts from the PSF. Details of the procedure will be given in a separate 
paper. We determined the PSFs down to 10~^ of the peak brightness both at the L15 
and L24 bands (Figs. [TJ) and e). The PSFs of the IRC MIR images are extended over the 
entire field-of-view due to the scattering within the detector array as recognized in the 



IRAC onboard Syi tzer ( IFazio et al. 



20041 ). which uses the same type of the detector array 



( IReach et al...2005. ). Correction for the artifacts was applied for each exposure frame after 



the sky background had been subtracted. Except for this additional process, the data were 
reduced by the standard imaging process with the latest flat data and the dithered images 
were co-added to obtain the final images (Figs. [IK and d). The obtained PSFs were fitted 
to the inner 20^' region of the co-added images and subtracted (Figs. [Tfc and f). The scaling 
of the PSFs is determined with an accuracy better than 0.5%. Extended structures due 
to the PSF patterns produced by the central star are cancelled out almost completely and 
the detached shell appears clearly in the central-star-subtracted images of both bands. The 
PSF well fits the image both at the inside and outside of the shell, confirming the detection 
of the sheU at L15 and L24. 



U Ant was also observed with the Far-Infrared Surveyor (FIS, 



Kawada et al. 



2007b 
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onboard AKARI on 2006 December 21 with the slow scan mode for a compact source. Data 
in four FIR bands, N60 (65 //m), WIDE-S (90 /xm), WIDE-L (140 /xm), and N160 (160 /xm) 
were simultaneously ob tained for a 10^ x 40^ a rea. The FIS data were processed in the same 



manner as described in 



Izumiura et al. 



( 120 lOl ) and the same PSF subtraction procedure was 



applied as for the IRC data. The FIS data are in agreement with the IRAS and PACS data 
within the uncertainties and have finer spectral sampling than them. The observational 
data used in this study, including the calibration data, are summarized in Table [H 



Results 



The central-star-subtracted image of L24 shows a smooth, symmetrical shell. The shell 
is also clearly seen at L15 (Figs. [U^ and f). The presence of an extended component is 
evident even in the original L24 image (Fig. [l]!), when compared to the PSF (Fig. [1^). The 
brightness peaks at around 41^' in the both bands. To investigate the shell structure in detail, 
we assume that the shell is optically thin at MIR and spherically symmetric. We integrate 
the fiux density in each projected annulus to obtain the 1-D radial brightness profile 
(RBP). The density di stribution of the shell along the radial direction is approximated by a 



Gaussian according to 



Maercker et al. 



( 120 101 ). The model emission is integrated along the 



line of sight, convolved with the PSF, and compared with the observed profile. We found 
that any single Gaussian distribution cannot give a reasonable fit with the observations 
and that the d o uble shell with the Gaussian parameters of shell 3 and 4 estimated in 



Maercker et al. 



( 120 lOl ) fits the observations fairly well. 



Figures [2^ and b show the double-shell model fit. The model fit parameters are 
summarized in Table [21 The errors in the RBP are estimated from the variance of the 
intensity in each annulus, which includes the error in the central star subtraction, the 
background sky fiuctuation, and the deviation from the spherical symmetry. Although the 
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MIR images do not resolve the second shell, the tail of the RBP at larger radii requires a 
second Gaussian component whose peak is located at around 50^'', especially at L24. The 
flux densities of both shells are estimated from the Gaussian parameters. The total shell 
fluxes at the FIS bands are estimated by the same procedure as for the IRC. The flux 
density of each shell cannot be determined separately at the FIS bands because of the low 
spatial resolution. 



Together with the FIS and PACS (IKerschbaum et al 



20 id ) data, we now obtain the 



spectral energy distribution (SED) of the detached shell of U Ant from 15 /xm to 160 //m, 
from which the temperature and mass of dust grains in the shell can be estimated. We 
assume that the variability at IR wavelengths is negligible because the visual magnitude 
amplitude is only about 1 mag. First we assume that dust grains in the both shells have 
the same temperature. The PACS data are given only for the sum of the flux densities of 
the central star and the shells, whereas the IRC and FIS data are available for the central 
source and the shells separately. Thus we estimate the dust temperature and mass by 
minimizing 

1=1 ^ 1=1 ^ 

^ (F^(A.) - A\f -C\fB^[\^T)f ^ 
1=1 

where B^{\^T) is the Planck function of the temperature T, F^(A^), F^(A^), and F^{Xi) are 
the observed flux densities of the central source, the shells, and the sum of the central source 
and the shells, respectively, at the wavelength i: i=l, 2, 3, 4, 5, 6, 7, and 8 correspond to 
15, 24 (IRC), 65, 90, 140, 160 (FIS), and 70, and 160 //m (PACS), respectively. The flux 
densities of the central source are estimated from the PSF flt, which include those from 
the star and the circumstellar envelope produced by the present mass-loss. We assume 
that flux densities are in the Rayleigh- Jeans regime in the spectral range in question. The 
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dust emissivity is assumed to be given by A~'^. In Eq. ([Tj) , T, A, and C are the fitting 
parameters, where T is the dust temperature and A and C are the scahng factor for the 
central source and the shell emission, respectively. In the fitting, we take account of the 
color corrections of each filter band (see Fig. [3D. 



The dust mass M can be estimated from 



M 



k(A) 



(2) 



where D is the distance to the star (260 pc) and is the dust mass absorption coefficient. 
The dust mass given by Eq. ([2]) is independent of the dust size when the dust size is much 
smaller than the wav elength in question. We adopt the dust parameters employed in 

( 119971 ) as /3 = 1.4 and /^(60 /xm)=150 cm^ g~^. They are i n agreeme nt with 



Izumiura et al. 



those of amorphous carbon grains with the specific density of 1.5 gem ^ (ISuh 



2000). 



The best fit result for the single-temperature case is shown in Fig. [3^, which gives 
T = 75.411^ K and M = 4:.3l:li x 1O"^M0. The single-temperature model cannot fit the 
observed SED well. The reduced is 6.6, which arises mostly from excess emission at 
15 jam. If the data point of 15 /xm is removed from the fit, the reduced improves to be 
1.9 and the model underestimates the observed fiux density at 15 /xm by a factor of 3. 

Alternatively, we assume that dust grains in eac h shell (shell 3 and 4) have different 



Maercker et al. 



f l2010h (Table Ej) are 



temperatures. The same shell parameters adopted by 
assumed. Using the fiux densities obtained for each shell at the IRC bands, we derive 
the best fit parameters for the two-shell model in a manner similar to Eq. ([T]). The best 
fit result for the two-shell model is shown in Fig. [Hb and its parameters are given in 
Table [21 The reduced becomes 1.5. The fit is significantly improved compared to the 
single-temperature model. The model fiux density at 15 /xm agrees with the observation 
within the uncertainty. The two-shell model well accounts for the observed SED from 15 to 
160 /xm. 
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Discussion 



Op tical scattered lig;ht observations indicate the presence of 4 detached sheUs around 



U Ant (IGonzalez Delgado et al. 



2001 



20031 : 



Maercker et al. 



20101 ) , among which the present 



observations detect MIR emission from shells 3 and 4. S hells 1 and 2 are seen only at 



the bands that contain the NaD lines and not full shells ( IMaercker et al. 



2010l ). It is not 



clear whether these are real shells or merely patchiness present in shells 3 and 4. Optical 
polarization images of dust-scattered light, line-scattered light data, and medium-resolution 
CO radio line observations suggest that shell 3 is gas-rich and most dust grains reside in 
shell 4. Our results indicate th at sh ell 4 dominates in the dust mass (Table [2j), confirming 



their results. 



Izumiura et al. 



( 119971 ) detect FIR emission from the detached shell at ^ 50^' 
at 60 and 100 /xm as well as that at ^ 3' around U Ant. They estimate the dust mass 
associated with the shell at ^ 50^'' as 2.2 x 10~^Mm , whic h is in good agreement with the 



present estimate of 1.6 x 10 ^M©. 



Maercker et al 



(I2OIOI ) derive the gas mass in shell 3 



as 2 X 10~^ Mq from CO observations and the dust mass in shell 4 as 5 x 10~^ Mq from 
polarization observations. The dust mass in shell 4 derived in the present study agrees with 
their result, taking account of the large uncertainty in the estimation from polarization 
observations. The dust-to-gas ratio in shell 3 is thus 8 x 10~^, supporting that shell 3 is 
dust-poor. The total dust-to-gas ratio for the sum of shells 3 and 4 is 0.008, which is in a 
reasonable range for carbon-rich AGB stars. 

Adopting the distance of 260 pc, the locations of shells 3 and 4 are 1.7 x 10^'' and 
1.9 X 10^'' cm from the star, respectively. Th e dynamical a^e of sh ell 3 is estimated as 2700 



yr from the CO line velocity of 19.5 km s ^ ( IMaercker et al. 



2010). The width of sheU 3 



indicates that the average mass-loss rate was 1.2 x 10 ''MQyr ^ at the shell formation. It 
should be noted that this includes the swept-up material and the shell may have widened. 



Detached shells around U Ant are resolved in PACS maps at 70 and 160 /xm. The 
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peak intensity is located at around 43'^ and they attribute it to shell 3 (IKerschbaum et al 



2010l ). The present analysis suggests that the FIR emission arises mostly from shell 4 
with a negligible contribution from shell 3. We calculate the RBP of the FIR emission by 
adopting the same radial dust density distributions as required to explain the MIR emission 
(Table Ej) and assigning the flux densities derived from the SED flt to each shell. Then it is 
convolved with the model PSF of PAC^. Figure [2]3 plots the result. The intensity peaks 
at around ^ 44^' and shows a broad proflle similar to the observ ations, which also suRges t 



2010|). 



the presence of an extended emission component beyond shell 3 ( IKerschbaum et al.l 
Because the column density peaks at the inner part of the shell on the line-of-sight and 
shell 4 has a relatively large FWHM, the peak of the RBP is shifted to a position inside the 
actual shell location. Therefore, the two-shell model explains the observed RBPs and SED 
from 15 to 160 /xm consistently. 



If amorphous carbon gr ains of 0.1 /im size 



ocated at 1.9 x 10^'' cm from the star are 



heated by a star of 5800 (jSchoier et al 



20051 ). their temperature will be about 70 K, 



which is marginally in agreement with the derived temperature of 60 K. On the other hand, 
the temperature of shell 3 (> 100 K) is deflnitely too hot since the estimated temperature 



at shell 3 is only 76 K. The dust dri 
two-s hell structure wa s suRR ested by 



't with respect to the Ra s as a n explanation for the 



on by 



Maercker et al. 



Gonzalez Delgado et al 



( 12003 1) and further elaborated 



( 120101 ). The estimated drift velocity is 3kms ^. In a steady state 
the drift velocity Vd is given by 



VgQpL 



cM 



(3) 



where Vg is the gas velocity, L is the stellar luminosity, M is the mass loss rate, c is the 
velocity of light, and Qp is the radiation pressur e efficiency factor, which is proportional 



to a, if a is much smaller than the wavelength (.Habing et al. 



19941 ). Thus small grains 



^http : //dirty . as . arizona . edu/'^kgordon/inips/conv.psf /conv_psf s . html 



have a low drift velocity. If grains of 0.1 fim have a drift velocity of 3kms~^, that of 10 nm 
grains will be lkms~^. This makes only a 2^' offset in the dust distribution relative to gas 
in 2700 yr, which cannot be resolved clearly in the present observations. Carbonaceous 
grains of 10 nm have small heat capacity and will be heated up to 100 K if they a bsorb a 
photon of leV by the temperature fluctuation mechanism (e.g., lAannestadl Il989l ) . Thus 
shell 3 can still retain very small grains, whose temperature can be quite high. These crude 
estimates suggest a possibility that smallest grains are still trapped in the gas-rich shell 3 
and emit hot MIR radiation. The actual relative motion between the gas and dust and the 
temperature of very small grains need to be investigated by numerical simulations for a 
given size distribution of dust grains. 



This work is based on observations with AKARI^ a J AX A project with the participation 
of ESA. The authors thank all the members of the AKARI project and AKARI stellar 
program members. This work utilizes the IRC data taken during the performance 
veriflcation phase. They are grateful to the IRC team members for their continuous 
encouragement and useful comments. 

Facilities: AKARI. 
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Table 1: List of the observational data. 



Observation ID 


AOT^ 


Date 




Object 




Note"^ 


1710071.1 


IRC02 a;L 


2006-12- 


■20 


U Ant 




target 


1710072.1 


FISOl 1.0;15;70 


2006-12- 


•21 


U Ant 




target 


1711365.1 


IRC02 a;L 


2007-05- 


•12 


IRAS 22396- 


■4708 


PSF & SCL 


1710039.1 


FISOl 0.5;15;70 


2006-11- 


•14 


IRAS 22396- 


■4708 


PSF 


4080019.1 


IRC03 a;L 


2007-03- 


•05 


CRL618 




SCL 


5124089.1 


IRC03 a;L 


2007-06- 


•19 


IRAS F06009- 


-6636 


PSF 


5124090.1 


IRC03 a;L 


2007-06- 


•22 


IRAS F06009- 


-6636 


PSF 


5124091.1 


IRC03 a;L 


2007-06- 


•22 


IRAS F06009- 


-6636 


PSF 


5124092.1 


IRC03 a;L 


2007-06- 


•21 


IRAS F06009- 


-6636 


PSF 


5124097.1 


IRC03 a;L 


2007-07- 


•01 


IRAS F06009- 


-6636 


PSF 


5124098.1 


IRC03 a;L 


2007-07- 


•04 


IRAS F06009- 


-6636 


PSF 


5124099.1 


IRC03 a;L 


2007-07- 


•04 


IRAS F06009- 


-6636 


PSF 


5124100.1 


IRC03 a;L 


2007-07- 


•05 


IRAS F06009- 


-6636 


PSF 


5124105.1 


IRC03 a;L 


2007-07- 


•12 


IRAS F06009- 


-6636 


PSF 



^Astronomical Observation Template for the IRC and FIS observations. See ASTRO-F Observer's Manual 
for details of the parameters (http: //www. ir . isas . jaxa. jp/AKARI/Observation/ObsMan/|). 

^PSF indicates the data that were used to derive the PSF. SCL indicates those used to estimate the scattered 
light contribution (artifacts). 
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Table 2: Parameters of the shell. 



Shell 3^ SheU 4^ Total 



Radial brightness profile (RBP) model 

Center position^ 43^/5 49^.7 — 

FHWM^ 2'!2 &!5 — 
Observed fiux density (Jy) 

IRC/L15 0.68 ±0.11 0.16lJ:^ 0.84 ± 0.05 

IRC/L24 2.30 ± 0.22 2.37 ± 0.23 4.67 ± 0.22 



FIS/N65 — — 25.8 



+5.3 
-3.1 



FIS/WIDE-S — — 20.1 ±4.2 

FIS/WIDE-L — — 8.4 ±3.1 

FIS/N160 — — 3.5 ±2.2 

PACS70 — — (24.21^:^)^ 

PACS160 — — {5'8toiY 
SED model results'^ 

Temperature (K) 103.6 ± 7.4 QO.OtH — 

Dust mass (M©) ^-^-0^72 x 10"^ ^-QOt^fg x 10"^ ^-Q^tofg x 10" 



^The designation of the shells follows 
^The RBP parameters are adopted from 



Gonzalez Delgado et al. 



(|200lh and 



Maercker et al. 



(|2010h . 
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^Estimated values by subtracting the stellar flux density from the observations (see 
^Results of the two-shell model (see 
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Fig. 1. — IRC MIR images of U Ant. The top row shows (a) the co- added image corrected for 
the artifacts, (b) the artifact-removed PSF pattern, and (c) the co-added image subtracted 
by the central star PSF for the L15 band. The bottom row (d, e, and f) shows the same 
images for the L24 band. The color scales are in units of ADU per pixel and the scale of the 
central-star-subtracted images is enhanced by a factor of 20 and 10 for the L15 (c) and L24 
(f), respectively. For c and f, the central 20'^ region is not plotted to show the shell structure 
clearly. 
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Fig. 2. — Radial brightness profiles (RBPs) of the detached sheU of U Ant at (a) 15 /im, (b) 
24 yum, where the profiles of the best-fit double-shell model are also indicated by the solid 
lines together with the RBPs of the two shells by the dotted and dot-dashed lines (see text 
for details), and (c) model profiles for PACS 70 (solid line) and 160 /im (dashed line) based 
on the two-shell model (see 




Wavelength (jim) 

Fig. 3. — SED fit results: (a) single-temperature and (b) two-shell case. The green solid 
line indicates the total flux density of the best-flt model and the red dashed and blue dot- 
dashed lines indicate the shell and central source flux densities of the model, respectively. 
The light blue dashed and purple dotted lines in (b) indicate the models of shells 3 and 
4, respectively. The black symbols with error bars indicate the observations, whereas the 
color symbols indicate the color-corrected model values of the lines of the same color, which 
should be compared with the observations. The circles and the squares indicate the shell 
emission (FJ), the diamonds the central source emission (F^), and the triangles the total 
flux densities (F^) (see Eq. ([I])). 



